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Abstract. We investigate the effect of supersonic relative velocities between baryons and dark matter, re- 
cently shown to arise generically at high redshift, on baryonic acoustic oscillation (BAO) measurements at 
low redshift. The amplitude of the relative velocity effect at low redshift is model-dependent, but can be pa- 
rameterized by using an unknown bias. We find that if unaccounted, the relative velocity effect can shift the 
BAO peak position and bias estimates of the dark energy equation-of-state due to its non-smooth, out-of -phase 
oscillation structure around the BAO scale. Fortunately, the relative velocity effect can be easily modeled in 
constraining cosmological parameters without substantially inflating the error budget. We also demonstrate 
that the presence of the relative velocity effect gives rise to a unique signature in the galaxy bispectrum, which 
can be utilized to isolate this effect. Future dark energy surveys can accurately measure the relative velocity 
effect and subtract it from the power spectrum analysis to constrain dark energy models with high precision. 
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1 Introduction 

Understanding the nature of dark energy is one of the key questions in contemporary science, and various 
techniques have been developed over the past decade to advance this goal. Measurements of the baryonic 
acoustic oscillation (BAO) feature in galaxy surveys provide an especially promising way to constrain the 
expansion history of the universe and the behavior of dark energy [1]. The physical scale of the acoustic 
oscillation is determined with high accuracy by cosmic microwave background (CMB) measurements (e.g., 
[2, 3]), making BAO measurements at low redshift a robust standard ruler Since current and future dark 
energy surveys aim to measure the BAO peak position at sub-percent level precision, significant efforts have 
been devoted to modeling nonlinear evolution and its effect on the BAO peak position [4-11]. At low redshift, 
where nonlinear effects are substantial on BAO scales, the BAO peak contrast is degraded and the peak position 
may be subtly shifted compared to the linear theory prediction. However, these nonlinear effects tend to 
produce smooth, broad-band power and may readily be removed by template fitting or explicit modeling of 
these effects [12-14]. 

Recently, Tseliakhovich and Hirata [15] discussed a new nonlinear effect in the growth of small scale 
structure at very early times. Prior to the recombination epoch, baryons are tightly coupled to the photons, while 
dark matter is decoupled from the baryon-photon fluid and its fluctuation grows logarithmically with small 
velocity. After cosmic recombination, when the tight coupling of the baryon-photon fluid is broken and the gas 
sound speed plummets, the relative velocity between the baryon and the dark matter fluids becomes supersonic. 
These random relative velocities effectively increase the Jeans mass for the formation of the earliest baryonic 
structure, thereby suppressing their abundance [ 15-17]. High-resolution numerical simulations [18, 19] have 
confirmed that supersonic relative velocities raise the minimum halo mass that can form the first stars, and 
delay star formation in time, by an amount that depends on the magnitude of the relative velocity (see also 
[20]). This modulation of baryonic structures at high redshift imprints signatures of the relative velocity effect 
in the power spectrum of objects such as minihalos at high redshift [15, 16], and hence the power spectrum of 
any observable which traces these objects can exhibit significant departures from simple linear biasing of the 
matter power spectrum on BAO (~ 100 h^^Mpc) scales. 

Even at low redshift z < 5, when massive dark matter halos form, the relative velocity effect can be 
indirectly important, contrary to the naive expectation that this effect becomes negligible within halos whose 
velocity dispersion greatly exceeds the relative velocity. One possible scenario [16] is that patchy reionization, 
driven by early minihalos, inherits the imprints of the relative velocity effect in the spatial distribution of early 
minihalos, resulting in a spatial variation of the subsequent star formation history modulated by the large-scale 
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power of the relative velocity effect. Consequently, massive galaxies at low redshift may retain the memory 
of the large-scale relative velocity effect seen in the minihalos at high redshift. Even tiny ^1% level effects 
on the colors or luminosities of massive galaxies can have significant implications, since these effects will be 
correlated across ~ length scales. Any residual relative velocity effect in low-redshift massive 

galaxies may have significant impact on the precision measurements of the BAO peak position. Since the power 
spectrum of relative velocities has a prominent oscillation structure that is out-of -phase with the matter power 
spectrum, it cannot be removed by a blind broad-band fitting of the galaxy power spectrum, and therefore it 
can potentially bias the determination of the BAO peak position. 

The primary purpose of the present work is to investigate the impact of any relative velocity effect that 
persists in low-redshift massive galaxies, on our ability to probe the underlying cosmology using galaxy power 
spectrum measurements, and to find possible ways to isolate and remove this effect in the analysis. The 
organization of this paper is as follows. In section 2 we briefly describe the supersonic relative velocity effect 
and its impact on the earliest baryonic structure. In section 3 we model the observed galaxy as a tracer of 
the underlying matter and the relative velocity distributions and compute the galaxy bispectrum to isolate 
the relative velocity contribution. We then investigate the impact of the relative velocity effect on the BAO 
measurements in section 4. We first compute the full power spectrum including the relative velocity effect in 
section 4.1 and quantify the shift in the BAO peak position due to the relative velocity effect in section 4.2. 
In section 5 we perform a Fisher matrix analysis to investigate the impact of the relative velocity effect on 
the cosmological parameter estimation and how well the relative velocity effect can be measured by using the 
bispectrum. Finally, we conclude in section 6 with a discussion of further implication. Technical details of our 
derivations are presented in Appendices A and B. 

Here we present our calculations assuming a flat ACDM universe with the matter density lj„i = ^mh^ = 
0.134 (17,„ = 0.271), the baryon density uj^ = Ofe/i^ = 0.0222 (f7f, = 0.045), the spectral index 7i,, = 0.966 
and its running = of the primordial curvature power spectrum = 2.42 x 10^^ (cts = 0.81) at 
Uq = 0.002Mpc^^ (consistent with the Wilkinson Microwave Anisotropy Probe measurements [2, 3] and the 
Sloan power spectrum measurements [21, 22]). 

2 Supersonic relative velocity effect 

Before the recombination epoch, baryons and photons are tightly coupled with sound speed Cg ~ l/VS, 
while dark matter is decoupled and cold. At the release of baryons following cosmic recombination, the gas 
sound speed drops precipitously (Cs ~ 10^^), rendering the baryon fluids supersonic (7M ~ 2 — 5). The 
different velocities of the dark matter and the baryon fluids result in suppression of the matter power spectrum 
at the Jeans scale (fcj ^ 200 Mpc^^) and suppression of the growth of structure at the relative velocity scale 
{kr = kj/Ai, A/h — IO^Mq ^ k~^) [15, 17]. Furthermore, the relative velocity effect modulated by large- 
scale acoustic oscillations imprints its signature in the collapsed baryon fraction at early time, and this effect 
can be used to probe the nature of minihalos before reionization [ 16]. 

The velocities ,„ of the baryon and dark matter distributions are related to their density fluctuations as 

k • k • 

Vb,,n(k,z) = «p a 5b^,„{k,z) = i— a Tb,,n{k,z) ipp{k) , (2.1) 

where ipp (k) is the primordial curvature perturbation in the total comoving gauge. The evolved density fluc- 
tuations are computed using the Einstein-Boltzmann code CMBFAST [23] and are expressed in terms of their 
transfer functions Jf, „j(k, z) = Tb.m{k, z) <Pp(k). Instead of taking the derivatives of the transfer functions, 
we simply use the velocity transfer function T„(fc, z) from the CMBFAST code (it is related to the density trans- 
fer function as Tyt,,vm = —o-T^.m/k). The relative velocity of the baryon and the dark matter fluids is then 
Vr = Vf, — v,„, and we define the dimensionless relative velocity as = Vr/cr,.„ (and their transfer func- 
tion Tru as well) using the one-dimensional rms relative velocity fluctuation cr,.„. Their statistical properties 
are completely determined by the two-point correlation function 

i^,,{T) = «(x)ui(x + r)) = i'^ir) 6,, + [V^,, (r) - ^x(r)] <5,. 6,. = / e"'"- P^f (k) , (2.2) 
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Figure 1. Transfer functions and correlation functions, (a) Relative velocity transfer function Tru{k) at 2; = 0. (fe) Ratio of 
the relative velocity to the dark matter transfer functions kTm/Tm (bottom panel). For comparison, the upper panel shows 
the ratio of the no-wiggle power spectrum to the linear matter power spectrum. The transfer functions are defined with 
respect to the primordial curvature perturbation, and the no-wiggle power spectrum is computed by using the Eisenstein 
and Hu [24] fit to the smooth power spectrum shape. 



where the indices i, j represent the spatial component of the relative velocity vector, the z-direction is set 
parallel to the separation vector r {'4>±_ = il}^! = ^yy, V^n = '^'zz), and the relative velocity power spectrum is 
defined as 

(4(kJ<(kb)) = (27r)3j^(ka -t- kfa) P:J(k„) (2.3) 



with 



P.ik) = ^ T^ik) P,{k) = ^ ^ . (2.4) 



For notational simplicity, we suppressed the time dependence of the transfer functions and their power spectra. 

Figure la shows the transfer function Tru{k) of the relative velocity effect at z = 0, where its amplitude 
indicates the linear growth of power given the primordial curvature power spectrum A^(fc) = A(fc/fco)"°~^ 
with A = 2.4 X 10^^ and fco = 0.002 Mpc^^. The relative velocity is predominantly sourced by structures 
on size scales of order the acoustic scale; larger-scale structures make little contribution, and the transfer 
function rru(fc) declines rapidly at fc < 0.03/iMpc^^ since baryons and dark matter behave similarly on scales 
beyond the sound horizon at decoupling. The ratio of the relative velocity to the dark matter transfer functions 
Tru{k)/Tm{k) is shown in Fig. lb, where the scaling factor 1/k reflecting energy-momentum conservation 
is removed. The prominent oscillation structure displayed in Fig. la is still visible in the bottom panel of 
Fig. lb, when compared to the matter transfer function. To facilitate the comparison of the acoustic oscillation 
amplitude and its phase in the relative velocity and the dark matter transfer functions, we compute the ratio 
of the no-wiggle power spectrum to the linear matter power spectrum and plot its square root in the upper 
panel of Fig. lb, where the no-wiggle power spectrum is a fit to the smooth shape of the linear matter power 
spectrum [24]. The oscillation structure in the relative velocity and the matter distributions is out-of-phase, 
and its amplitude (~5%) is smaller in the matter distribution as the acoustic oscillation was only present in the 
baryon distribution, not in the dark matter distribution at early time. By contrast, the relative velocity arises 
entirely from the acoustic oscillations, and hence its oscillation amplitude is fractionally order one. 

However, the relative velocity of the baryon and the dark matter distributions is not directly observable. 
Dalai et al. [16] argue that the collapsed baryon fraction at early time is affected by the relative velocity effect 
as it changes the effective sound speed and increases the equivalent Jeans mass for baryonic gas to collapse in 
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dark matter halos. The two-point correlation function of the collapsed baryon fraction can be computed [16] as 

cx (u^(x)w,2,(x + r)) , (2.5) 

implying that the collapsed baryon fraction is a biased tracer of the relative velocity w^(x). We quantify the 
effect of the relative velocity on the BAO peak shift in section 4. 



3 Separating relative velocity contribution: Bispectrum 

At high redshift, z > 10, the supersonic relative velocity allows baryons to advect out of small dark matter 
halos, effectively increasing the Jeans mass of the gas in a velocity-dependent manner, thereby modulating 
the large-scale clustering of collapsed baryonic objects [ 16]. Additionally, the abundance of dark matter halos 
is also modulated by a similar effect [15]. Quite generally, we can write the large-scale fluctuation of the 
collapsed objects as [16] 

(5g(x) bi (5m(x) + br [w^(x) - a^„] , (3.1) 

where bi and br are the bias parameters of the collapsed halo with respect to the matter density and the relative 
velocity, respectively. Hereafter, we collectively call the collapsed baryon and dark matter system a galaxy, 
though a galaxy is often used to refer to a baryon only system and at high redshift the collapsed baryon and 
dark matter systems may appear different from typical "galaxies" at low redshift. The galaxy fluctuation is 
not only a tracer of the underlying matter distribution, but also a tracer of the relative velocity. The magnitude 
of the relative velocity bias br is computed by Dalai et al. [16] at z > 10, where the collapse of the baryonic 
structure can be relatively simply modeled. As discussed above, the relative velocity between baryons and dark 
matter can indirectly modulate the properties of galaxies at low redshift. The large-scale clustering properties 
of observed galaxies can, to lowest order, therefore be written in the form of eq. (3.1). In this case the velocity 
bias br is treated as a free parameter, reflecting our great uncertainty in the physical processes that determine 
galaxy properties at low redshift. 

Since the relative velocity Ur enters quadratically in eq. (3.1), it is readily apparent that the relative 
velocity effect will lead to a nonvanishing bispectrum of the galaxy fluctuation, providing a direct way to 
isolate its contribution from the matter density distribution. Noting that the relative velocity has no directional 
dependence, the ensemble average of the galaxy fluctuations can be computed as 

/^3q 
7^~T3 X! ('5g(ka)(5g(k,,X(q)'u*.(kc - q)) + cyclic permutations (3.2) 

i—x,y,z 



I J7^ I J7^ I 7;^5I("^(^i)^^(^''~^i)^^(^2X(k6-q2)w!!(q3)u!?(kc-q3)) 



. (27r)3 J (2^)3 J (2^)3 

Counting only the connected part (k^ + kb + kc = 0) in eq. (3.2), we can derive the galaxy bispectrum (see 
Appendix A) 

B,{K,kM^2blbr J2 [n(ka)n(M+cycl.] +8 6^ J ^^/>^^(k,+q)Ff (k,-q)i-„^'(q) 

i—x.y,z i-j^^ 

)] (3.3) 

/(Pa 
J^PrnilK + q|)-Pm(|kb - q|)P,„ (<?)G„ (k^ + q, kfc - q)G„(ka + q, q)G„(q, kfe - q) , 



d3q 

where the cross-power spectrum of the matter density and the relative velocity is defined by 

(5™(k„)<(kf,)) = (2^)35^(k, +kb) PJ(k„) = -K(ka)5™(kb)) (3.4) 

as 

Pi(k) T„,{k)P^{k) = ^ ^ Tru{k)Trn{k)P^{k) = ^ ^ ^^P,n{k) , (3.5) 
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Figure 2. Relative velocity contributions to the galaxy bispectrum for an equilateral (left) and a squeezed isosceles (right) 
triangular shapes. Dotted, dashed, and dot-dashed lines show each component of the galaxy bispectrum in eq. (3.3), and 
their sum is shown as solid lines. The bias parameter br/bi = 0.01 is assumed for computing the cubic term in eq. (3.3). 



and the relative velocity kernel is 



G (k k ) T^vika) Ty^ki,) k^ ■ k^ T^uika) Tru{kb) ^ 

with fXab = ka • kfc/fca/Cfc. 

Figure 2 illustrates the galaxy bispectrum for two triangular configurations. The left panel examines the 
scale-dependence of the galaxy bispectrum for an equilateral shape (ka ^ kf, ^ kc = k, /.i = —0.5). With 
the relative velocity kernel Gu{k, k) oc on scales of interest, the galaxy bispectrum declines sharply on 
small scales in proportion to k~^-^, while it flattens on large scales. However, because baryons and dark matter 
have similar velocities on large scales k < seen in Fig. la, this term in the galaxy bispectrum 

asymptotically vanishes. The acoustic oscillation structure of the galaxy bispectrum around k ~ 
reflects the oscillation due to the relative velocity effect. Being a convolution, the cubic contribution (dot- 
dashed) in eq. (3.3) is constant and nonvanishing on large scales, because the contribution to the cubic term 
arises at g ~ 0.01 — 0.1/iMpc~^. The right panel illustrates the scale-dependence of the galaxy bispectrum for 
a squeezed isosceles triangle [ka = fcf, = fc 7^ fee, fiab — —0.99). Since the third scale of the isosceles is larger 
kc = k^2 + 2yU = 0.14fcandthe scaling of the bispectrum is steep, two contributions P,„(fc)P„i(A:c)G„(fc, fee) 
in eq. (3.3) that mix different scales k and fee are larger than the other term P„i {k)Pjn{k)Gu{k, k) on small 
scales, but it falls over faster on large scales. The cubic term as in the equilateral configuration is subdominant 
fovbr/bi < 0.1 at A: > 0.01/iMpc"\ 

At low redshift, however, the matter fluctuation develops substantial nonlinearity on small scales and 
is quasi-linear even at relatively large scale k ~ 0.1/iMpc~^, demanding treatment beyond hnear order in 
eq. (3.1). The nonlinear evolution in the matter density distribution results in a nonvanishing bispectrum of 
the galaxy fluctuation even in the absence of the relative velocity effect. Here we adopt standard perturbation 
theory to compute the nonlinear terms in the matter density distribution to the third order (e.g., [25-27]): 
dm (k) = (5m ^ (k) + Sm^ (k) + (5m ^ (k), where the superscript indicates the perturbation order and the second- 
order matter fluctuation is 

€'(k)= I ^<5W(q)<5W(k-q)F2(q,k-q) (3.7) 
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Figure 3. Scale-dependence of the full galaxy bispectrum for two triangular shapes. Three different components contribute 
to the full galaxy bispectrum in eq. (3.10); The nonlinear evolution of the matter density distribution (dot-dashed), the 
nonlinear bias (dotted), and the relative velocity effect (dashed). The cubic term in eq. (3.10) is omitted to avoid clutter. 
The full galaxy bispectrum is shown as solid lines. The bias parameters 62/61 = 0.1 and br/bi — 0.01 are assumed. 



5 2 / ka • k;, \ ka • kf, / 1 



with the second-order kernel 

Along with the nonlinear gravitational evolution in the matter density distribution, we also need to con- 
sider nonlinear galaxy bias. Expanding to third order, we model the galaxy fluctuation as 



(x) = 61 S,n (x) + ^ [51 (x) - al] + ^ 51 (x) + K (x) - J , (3.9) 



where we keep only the linear order term in the relative velocity effect as it is the relic effect of the early 
universe and it decays with the expansion factor a, rendering higher-order terms in w^(x) negligible at low 
redshift. Therefore, the full bispectrum of the galaxy fluctuation in eq. (3.9) is 

Sg(ka,kb,ke) = bl [2P^{ka)P,n{h)F2{ka,kb) + Cycl] + ^ bj [P™ (fca)P™ (^b) + Cycl.] (3.10) 
+blbr [2P^{ka)Pm{h)Gu{kaM) + Cycl.] 

+8 bl J ^^^™(|ka + q|)Pm(|kb - q\)P„^iq)Guika + q, kb - q)G„(ka + q, q)G„(q, kb - q) , 

where Pm (k) is the linear matter power spectrum. In addition to the relative velocity contributions in eq. (3.3), 
the full galaxy bispectrum receives two more contributions from the nonlinear evolution: The first square- 
bracket represents the contributions of the nonlinear evolution in the matter density distribution, while the 
second square-bracket represents the nonlinear bias contributions described in eq. (3.9). 

Figure 3 plots the scale-dependence of the full galaxy bispectrum in eq. (3.10). The left panel dissects 
each component of the galaxy bispectrum given an equilateral triangular shape {ka ^ kb ^ kc = k, fi = —0.5, 
F2 = 0.29). In this and subsequent figures, we assume the nonlinear bias parameter 62/^1 = 0.1 and the 
relative velocity bias parameter br/bi — 0.01 as our fiducial bias parameters for illustration. Two contributions 
from the nonlinear bias (dotted) and matter density (dot-dashed) are comparable in amplitude if 62 — bi, and 
they are dominant over the relative velocity effect (dashed) on small scales. However, the relative velocity 
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Figure 4. Full galaxy power spectrum with the bias parameters &2/fei ~ 0.1 and br/bi — 0.01. The galaxy power 
spectrum (thick solid) is computed by using eq. (4.2). Thin solid, long dashed, and dot-dashed lines represent the auto- 
power spectrum of the relative velocity effect (~ b^), the cross-power spectrum of the relative velocity and the nonlinear 
galaxy bias (~ br), and the auto-power spectrum of the nonlinear galaxy bias, respectively. The latter (long dashed) 
becomes negative at fc = 0.04/iMpc^^, beyond which its absolute value is plotted. The linear and the nonlinear matter 
power spectra are shown as dotted and short dashed lines (the short dashed line is largely obscured by the thick solid line). 



contribution is comparable to the matter density contribution at large scale k ~ 0.03/iMpc~^, corresponding 
to the peak seen in Fig. la. The right panel shows the galaxy bispectrum for a squeezed isosceles shape 
{ka = kf, = k kc, f-iab = —0.99). For the squeezed configuration, the large scale power Pm{kc) with kc = 
0.14fc enhances all three components on small scales and reduces the components on large scales. The relative 
velocity peak leaves a bump at fc ~ in the squeezed configuration, similar to the scale seen in 

the equilateral case. Therefore, the unique signature of the relative velocity effect in the galaxy bispectrum on 
large scales can be used to robustly measure the relative velocity effect for br > 0.01 in low-redshift massive 
galaxies. 



4 Effect on BAO measurements: Power spectrum 

In addition to generating a large-scale bispectrum, the relative velocity effect can also modify the galaxy power 
spectrum on scales of order the sound horizon. Numerous ongoing and future galaxy surveys seek to measure 
the galaxy power spectrum on these scales with extraordinary precision to localize the BAO feature and thereby 
reconstruct the expansion history of the universe, constraining the kinematic properties of dark energy. From 
the viewpoint of BAO probes of dark energy, the relative velocity effect could be a significant contaminant, as 
we illustrate in this section. 



4.1 Relative velocity contribution and galaxy power spectrum 

We compute the power spectrum of the galaxy fluctuation using eq. (3.9), but for simplicity we first compute 
the relative velocity contributions only. With caution that the cross-power spectrum (fc) in eq. (3.5) depends 
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on the direction of the wavevector, the ensemble average of the galaxy fluctuation is 



(27r)3 J (27r) 



12_ \ - 

^3 



2 6ii^2(qi,ka - qi) ((5.m(qi)(5.m(ka - qiX(q2X(k6 - q)) 



+b2 {S,n{qi)6,n{ka - qiX(q2X(kb - q)) 



(27r)3 J (27r)3 



X! «(qiX(ka - qi)u;^(q2)-u^(k;, - q)) , 



i,J=x,y,z 

and by isolating the connected part only the galaxy power spectrum can be computed as 



P,(k) = kJ -0 ^PJ(q)Pi(k-q) 

i—x,y,z 



d^ 

( 

(27r)3 



4 6iF2(q,k-q) + 2 62 



i,j—x,y,z 



P„(g)P„(|k-q|)G4q,k-q) 



4 6i&.P2(q, k - q) + 2 foa^r + 2 fe^G„(q, k - q) 



(4.1) 



(see Appendix A for derivation). The third term in the square-bracket is the pure relative velocity contribution 
to the galaxy power spectrum, while the other terms represent the contributions of the relative velocity effect 
in conjunction with the nonlinear galaxy bias. 

Finally, accounting for all the remaining contributions from the nonlinear evolution in the matter density 
distribution, the full power spectrum of the galaxy fluctuation can be written as 



i&^ + 2&i62P2(q,k-q) 



P,(k) = 6?PNL(k)+ J Aq_p„^(q)p„(|k-q| 

+4 6i6,P2(q, k - q)G„(q, k - q) + 2 62&.G„(q, k - q) + 2 fe^G^lq, k - q)^ 



(4.2) 



To second order in the power spectrum, the linear bias parameter bi is renormalized. Pnl (fc) in the first term 
is the nonlinear matter power spectrum, while Pm{k) in the integral is computed by using the linear matter 
power spectrum [28, 29]. Compared to eq. (4.1), two additional contributions in the square-bracket arise from 
the nonlinear galaxy bias. 

Figure 4 shows the full galaxy power spectrum computed by using eq. (4.2). Assuming our fiducial bias 
parameters &2/&1 — 0.1 and br/bi = 0.01, the galaxy power spectrum (thick solid) is largely determined 
by the matter power spectrum (linear: dotted, nonUnear: short dashed) on all scales, yet the relative velocity 
effect and the nonlinear galaxy bias affect the galaxy power spectrum at the percent level on various scales. 
These contributions expressed in the second square-bracket of eq. (4.2) are split into three components with 
different dependences on the relative velocity bias br'- The auto-power spectrum of the relative velocity effect 
(b^: thin solid), the cross-power spectrum of the relative velocity and the nonlinear galaxy bias {br'. long 
dashed), and the auto-power spectrum of the nonlinear galaxy bias (dot-dashed). The auto-power spectrum 
(thin solid) of the relative velocity effect closely resembles the velocity power spectrum with the prominent 
oscillation structure seen in Fig. la. The auto-power spectrum (dot-dashed) of the nonlinear galaxy bias is 
constant on large scales and approaches the shape of the matter power spectrum on small scales, while its 
power is enhanced at fc ~ 0.2/iMpc^^ by the nonlinear kernel P2(q, k — q) that puts more weight on the 
large-scale power in the convolution. The cross-power spectrum (long dashed) of the relative velocity and the 
nonlinear galaxy bias changes sign at k ~ 0.04/iMpc^^, reflecting the out-of-phase nature of the velocity and 
the matter distributions seen in Fig. lb. 

Future galaxy surveys that use BAO probes of dark energy will measure the large-scale galaxy power 
spectrum and attempt to determine the BAO peak position at sub-percent level precision. At fc ~ 0.1/iMpc^^, 
the auto-power spectrum of the relative velocity effect (thin solid) is rather sub-dominant compared to the con- 
tribution of the cross term and the nonlinear galaxy bias term (long dashed and dot-dashed) for our fiducial 
parameters. Therefore, both the nonlinear bias 62 and the relative velocity bias b,. parameters must be deter- 
mined well in order to completely clean out the relative velocity effect in the galaxy power spectrum. This 
trend remains unchanged for br/bi < 0.1 given 62/&1 = 0.1. 
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Figure 5. Impact of the relative velocity effect on the BAO peak shift. Various lines represent the peak shift (a — 1) in 
percentage as a function of the relative velocity bias parameter br/bi, given the value of the nonlinear galaxy bias parameter 
h2/b\. Following the Seo et al. [12] approach, the peak shift is obtained by fitting the template power spectrum to the full 
galaxy power spectrum in eq. (4.2) accounting for the nonlinear growth and the anomalous power (see text for detail). 

4.2 Shift in the BAO peak position 

We quantify the impact of the relative velocity effect on the precision measurement of the baryonic acoustic 
oscillation feature in future dark energy surveys, assuming that no information is utilized for cleaning out the 
relative velocity effect in the power spectrum measurement. Following the Seo et al. [12] approach to modeling 
the scale-dependent nonlinear growth and anomalous power in the matter power spectrum, we fit the template 
power spectrum to the full galaxy power spectrum computed in eq. (4.2) to characterize the BAO peak shift. 
The template power spectrum is parametrized as 



with polynomial coefficients Ci and ai taken as free parameters. Any deviation of the parameter a from unity 
represents the BAO peak shift in measurements. The two groups of polynomials in k (multiplicative and addi- 
tive) account for the scale-dependent nonlinear growth and the additive broad-band power; This parametriza- 
tion corresponds to the "Poly7" fit in Seo et al. [12]. The evolved linear matter power spectrum is then given 
in the form 



where Pno- wiggle (fc) is the fit to the smooth power spectrum shape without the BAO wiggle from Eisenstein 
and Hu [24] and = 8.8 ft,~^Mpc accounts for the degradation of the BAO wiggle in time {z = 0) [12]. 

In Fig. 5 we perform a analysis to compute the BAO peak shift as a function of the bias parameters 
hr/hi and 62/&1 by fitting the template power spectrum in eq. (4.3) to the full galaxy power spectrum in 
eq. (4.2) over a range of < k < Since nonlinear effects do shift the BAO peak 

position (~ 0.5% at 2 = 0.3 [12]) at low redshift, we isolate the effect of the relative velocity on the BAO 
peak shift by setting Aa = when hr/bi = 0. The relative velocity effect shifts the BAO peak position no 
more than 1% at hr/bi < 0.01 with a reasonable range of the nonlinear galaxy bias parameter, but its impact 
increases dramatically at br/bi S> 0.01, because the acoustic structure in the relative velocity effect is anti- 
correlated with the structure in the matter distribution. While a large nonlinear galaxy bias can affect the BAO 



Ptik) = (cq + cik + C2fc^)Pcvo(fc/a) + (flo + ciik + asfc^ H h ayfe^) 



(4.3) 




(4.4) 
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peak position, its parameters are, by contrast, relatively well constrained, and its impact is at the sub-percent 
level for our fiducial value hr/hi over a range of 62/^1 values considered here. For negative values of h^jhx, 
the cross-power spectrum (long dashed in Fig. 4) of the nonlinear galaxy bias becomes positive, but the auto- 
power spectrum (dot-dashed in Fig. 4) of the nonlinear galaxy bias is reduced due to the sign change in 62- At 
k > the latter is dominant over the cross-power spectrum, and hence the peak shift is reduced 

at br/bi = 0.01, compared to the case with > 0. However, as the relative velocity bias increases, the 

cross-power spectrum contributes more to the peak shift, changing its direction. The auto-power spectrum (thin 
solid in Fig. 4) of the pure relative velocity contribution around the BAO scale is rather sub-dominant over the 
scales for the range of 6r/&i considered here. 



5 Cosmological constraining power and bias in parameter estimation 

In section 4 we showed that the presence of the relative velocity effect can adversely affect our ability to 
measure the BAO peak position and constrain the cosmological parameters when its effect is ignored in the 
power spectrum analysis. However, it was also demonstrated in section 3 that measurement of the bispec- 
trum provides a promising way to measure the relative velocity effect and account for it when we estimate 
the cosmological parameters. In this section, we perform a Fisher matrix analysis to answer three ques- 
tions; What is the bias in our parameter estimation if the relative velocity effect is unaccounted for? How 
well can we measure or constrain the relative velocity effect using the galaxy power spectrum and bispec- 
trum measurements? And how much would the cosmological parameter constraints be inflated if we included 
the relative velocity bias as additional parameter? To answer these questions and for definiteness, we con- 
sider a galaxy survey of volume V ^ 10 {h~^Gpc)^ with number density Ug = 10~"'(/7,~^Mpc)^. As our 
cosmological model, we assume a flat universe with a constant dark energy equation-of-state and add two 
bias parameters &2/&1 and br/bi to describe the nonhnear galaxy bias and the relative velocity distribution: 
p = (tt-s, as, o;;,, Wdc, wq, &2/&1, br/bi) with /i^ = + Wdc (see section 1 for our choice of the fiducial 
model parameters). ' 

First, we compute the impact of the relative velocity effect on cosmological parameter estimation from 
galaxy power spectrum measurements when the relative velocity effect is unaccounted. This case will represent 
the currently planned galaxy surveys, as the null hypothesis is devoid of the relative velocity effect. Here 
we quantify its impact in terms of the parameter bias Awo in the dark energy equation-of-state. Our model 
parameters are described by p = p \br/bi = 0], in which the relative velocity bias parameter is incorrectly set 
zero {bj. = 0), while in reality it is nonzero {br 7^ 0). The parameter bias in our estimation of the dark energy 
equation-of-state is then (see Appendix B) 



F-i(p) 



E 



1 dPgjk) 

i (fc) dp. 



Pg(fc|p)-Pg(fc|p) 



(5.1) 



where the Fisher information matrix is 



fcmax 

E 

fc=*;mi, 



dPg{k\f>) dPg{k\P) 



(fc) dp. 



dpj 



(5.2) 



the uncertainty in the power spectrum measurements is [30] 



Air^k^Ak F/(27r)3 



P,(fc) + - 



(5.3) 



and k„ 



Ak = 2Tr/V 



1/3 



0.0029 /iMpc . In computing the Fisher information matrix in eq. (5.2) we 



also add the Planck prior on the cosmological parameters, following the procedure described in the Appendix 
of the Dark Energy Task Force final report [I] (see also [31, 32]). 



'We assume that the hnear bias bi and the matter fluctuation normalization is degenerate in the power spectrum analysis, and we 
can only constrain its combination. The degeneracy is partially broken by the nonlinear effect and also by the bispectrum, but since the 
normalization is nuisance in the present analysis we combine it with the linear bias term and remove it from our parameter set by assuming 
that other parameters are not affected by the change in the overall normalization of the measurements. 
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Figure 6. Impact of the relative velocity effect on the cosmological parameter estimation. For computing the constraint 
on the dark energy equation-of-state wq, we use the galaxy power spectrum measurements in a survey of volume V — 
10 (/i~^Gpc)^ with number density Ug — 10~'*(ft~^Mpc)^'^ Parameter constraints are obtained by adding CMB priors 
from Planck and marginalizing over the remaining parameters (see text for details), (a) Parameter bias in measuring too 
when the relative velocity effect is unaccounted, (b) Measurement uncertainty in the relative velocity bias parameter (the 
nonlinear galaxy bias parameter is fixed = 0.1). (c) Fractional increase in the cosmological parameter constraints 

when the relative velocity bias is modeled and marginalized over (fcmax = 0.2/iMpc~^). 



Figure 6a shows the bias in the dark energy equation-of-state wq as a function of the relative velocity bias 
given various values of the nonlinear galaxy bias 62/^1- We have fixed fc,„ax — 0.2ft,Mpc^^. Naturally, the bias 
Awq increases with larger value of K/bi, as the relative velocity effect is unaccounted for Negative values 
of reduce the auto-power spectrum (dot-dashed in Fig. 4) of the nonlinear galaxy bias and change the 
sign of its cross-power spectrum (long dashed in Fig. 4) with the relative velocity on small scales, enhancing 
the difference in the power spectrum around the BAO scale and thereby increasing the bias Awq, compared 
to reducing the relative velocity contribution in the fiducial model with positive 62/61. Meanwhile, since our 
parameter estimation is affected by the power spectrum measurements on all scales rather than just on BAO 
scales, there is a cross-over point in br/bi, at which the bias from the cross-power spectrum is cancelled by 
the bias from the auto-power spectrum of the relative velocity. Moreover, the parameter bias obtained here 
is rather insensitive to the precision of the power spectrum measurements, since it arises from our incorrect 
modeling of the mean value of the power spectrum. For the fiducial value 62/61 — 0.1 (solid), the bias Awq 
is at the 10% level for 6^/61 = 0.01 and 2% for 6^/61 = 0.001. While a more thorough analysis of the BAO 
peak position may further reduce the parameter bias than our simple estimates using the Fisher matrix (like 
in section 4.2, e.g., [12, 33]), Fig. 6a demonstrates the significance of the relative velocity effect on the future 
dark energy surveys, and the relative velocity bias parameter should be included and marginalized over in the 
power spectrum analysis. 

Next, we estimate our ability to measure the relative velocity effect from the galaxy power spectrum 
and bispectrum measurements. Our fiducial model in this case is then described by p[6r/6i ~ 0], correctly 
representing the reality. The measurement uncertainty in the relative velocity bias parameter can be expressed 
as 

</.. = [^-'(P)] . (5-4) 



and the Fisher matrix Fij now includes that of the bispectrum when the bispectrum measurements are used 
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[34-36], 



E E E 



1 dBg{ka,h,kc) dBg{ka,kb,kc 



with the bispectrum variance 



'B. 



^abc 



8n^kakbk,Ak3 V/{27r)3 



dpi 



PM + - 

I if-, 



Pg{kb)+— 

n„ 



Pg{kc) + — 



(5.5) 



(5.6) 



where k*-^^^ = max(fcmin, \ka — kb\) and the symmetry factor is Sabc ~ 6, 2, 1 for equilateral, isosceles, and 
general triangular configurations, respectively. We have assumed that the bispectrum estimates are Gaussian 
distributed. 

Figure 6b shows the measurement uncertainty ObrlW ™^ 'he relative velocity bias parameter, after all 
the remaining parameters p are marginalized over, which can be translated into the detection significance 
SjN ~ {br/bi)/<J{b,./bi) ■ With the few percent level contribution of the relative velocity effect at 6,. /6i = 0.01 
around the BAO scale seen in Fig. 4, power spectrum measurements (solid) alone can constrain the relative 
velocity bias as small as br/bi = 0.01 at the 2 — 3(7 confidence level. Bispectrum measurements (dashed) 
provide additional but less stringent constraint on the relative velocity bias parameter if h,. = 0, than power 
spectrum measurements, especially when the basic cosmological parameters are already well constrained by 
CMB measurements. Hence, the combination (dotted) of the power spectrum and bispectrum measurements 
yields constraints mainly derived from the power spectrum measurements. However, as illustrated in Fig. 3, 
the bispectrum exhibits a unique signature of the relative velocity effect on large scales, if the relative velocity 
effect is present. Therefore, if any hint of the relative velocity effect is found in the power spectrum anal- 
ysis, measurements of the bispectrum can be used to confirm and constrain the relative velocity effect in an 
indisputable way. 

Finally, we investigate the cost in the parameter constraints by adding the additional parameter br/bi and 
marginalizing over it. Figure 6c compares the parameter constraints in two cases, when the relative velocity bias 
parameter is modeled, and when it is simply assumed zero, i.e., ap^ (p)/o'p„ (p) ■ The matter and the baryon 
density constraints auj^ and a^^ are little affected by the addition of the relative velocity bias parameter, as they 
are largely independent. The constraints on the spectral index Ug and its running as become weaker by a few 
percent level, compared to the constraints in the model without the relative velocity bias parameter, indicating 
that the relative velocity effect is somewhat degenerate with Ug and as- Last, as we showed in previous sections 
the relative velocity effect can adversely affect the dark energy parameter estimation, and Fig. 6c shows that 
the constraints on dark energy parameters wq and cjdo would degrade by ^ 8%, if we included the relative 
velocity bias parameter While the spectral index and its running are also affected by the relative velocity 
parameter, they are highly constrained by CMB measurements, in contrast to the dark energy parameters. 
Meanwhile, the inflated constraints on the dark energy parameters can be reduced by using the bispectrum 
measurement to improve the constraint on the relative velocity bias parameter and thereby recovering the dark 
energy constraints. Overall, the cost of modeling the relative velocity bias is rather low, especially considering 
the magnitude of the parameter bias when the relative velocity effect is unaccounted for 



6 Discussion 

We have investigated the relative velocity effect of the baryon and the dark matter distributions on the bary- 
onic acoustic oscillation (BAO) measurements at low redshift. The relative velocity effect imprinted in the 
earliest baryonic structure at high redshift leaves its signature in the power spectrum of these objects [16], 
which in turn alters the subsequent star formation and reionization history, affecting the large-scale clustering 
of massive galaxies at low redshift. The relative velocity effect that persists until the late time results in a 
nonvanishing galaxy bispectrum, and its contribution to the bispectrum can be easily isolated on large scales 
k ~ 0.03/iMpc~^, where the relative velocity contribution peaks. Since the acoustic structure of the relative 
velocity effect is non-smooth and out-of-phase compared to the structure in the matter distribution, its effect 
cannot be removed by using a broad-band filter that is designed to deal with the nonlinear matter evolution, 
and the relative velocity effect can potentially bias the BAO peak position determination. 



-12- 



We find that the relative velocity effect can, if unaccounted, shift the BAO peak position measurements 
and bias the estimates of the dark energy equation-of-state by ~ 10%, depending on the amplitude of the relative 
velocity effect. However, the power spectrum measurement itself is sensitive enough to the relative velocity 
effect that we can measure its amplitudes as small as hr/hi ~ 0.01 along with other cosmological parameters. 
Including the relative velocity bias as additional parameter may inflate the constraint (7^,^ on the dark energy 
equation-of-state by 8%, a marginal cost to pay, compared to the potentially pernicious risk arising from the 
unaccounted relative velocity effect. Bispectrum measurements also provide tight constraints on the relative 
velocity effect, though its constraint in the fiducial model without the relative velocity effect is less stringent 
than from power spectrum measurements due to larger measurement uncertainties. However, the bispectrum 
signal is substantially enhanced, especially on large scales, if the relative velocity effect is present, providing 
alternative way to confirm and measure the relative velocity effect in a more robust and model-independent 
way. 

Many future dark energy surveys such as BigBoss,^ Euclid,^^ and the Wide-Field Infrared Survey Tele- 
scope"^ plan on measuring galaxies at moderately high redshift, ^ > 1, where nonlinear growth is less severe 
and even smaller scales may safely be included in the analysis. However, these high redshift surveys are not 
immune to the relative velocity effect, which grows in amplitude like (1 + z), eventually dominating over 
the matter fluctuations at z 10. At the redshifts probed by these BAO surveys, the largest uncertainty in 
our calculation is the amplitude of the relative velocity bias parameter hr- Any effect of relative velocities on 
massive galaxies at such redshifts must be indirect, making difficult any theoretical estimation of the relative 
velocity bias parameter It is quite plausible that the complicated process of massive galaxy formation obliter- 
ates the relative velocity effect at low redshift, meaning that this effect may not be apparent in the distribution 
of low-redshift massive galaxies. Nevertheless, we note that the relative velocity effect can be easily modeled 
without significantly inflating the cosmological parameter constraints, and if present but unaccounted, it can 
catastrophically affect the ability of BAO survey to perform high precision cosmological parameter estimation. 
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A Power spectrum and bispectrum computation 



Here we derive key equations for computing the full galaxy power spectrum in eq. (4.2) and bispectrum in 
eq. (3.10). Full computation of the galaxy power spectrum requires a volume integration over the wavevector q. 
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and it can be obtained by a combination of the two integrations: 



/ (0j^™(9)^™(|k-q|)^^2(q,k-q) 



27r2 



3r + 7fi- lOr/i^ 



14 r(l + H - 2r/i) 



/ (^^™(9)^™(|k-q|)Gn(q,k-q) 



(A.l) 



(A.2) 



27r2 



Tru{q) Tru{ky/l + - 2r^) r - 



Tm(g) r„,(A:v/l + r2 -2r^) ^1 + r2 - 2r/i 



where r ~ q/k is the ratio of the two wavevectors k and q and /i is their cosine angle. 

For the bispectrum computation, since the connected wavevectors (ka + kf, + kc = 0) defines a 2D 
surface rather than a Une (ka + kb — 0) in the power spectrum case, the azimuthal symmetry in the volume 
integration is broken, and a full 3D integration needs to be performed. Given a triangular configuration {ka, kf,, 
and ^ab = ka • 'k.b/kakb), the last term of the full bispectrum in eq. (3.10) is 

J (27r)3'^™^''^'' ^ 'i\)Pr,ii\kb - q|)Pm(g)G„(ka + q,kb - q)G'„(ka + q,q)G„(q, k^ - q) 

r2jfcr) dfi 



(A.3) 



T,l{ka^l + r^+2rfi) 



^/l + 7-2 + 2r^ 



r.l(fcaV7-2+r^2_2rr->,,) 
T2 (^^^7^2 + ^,2 _2rrVfc,) 



i/l + r2 + 2ryL^r^ + r'2 — 2rr' \/r'^ + r'"^ — 2rr' 



with = fcb/ fed and /x^g = MabM + ^^1 — Mab\/l ~ cos , where is the azimuthal angle of q. 



B Parameter forecast and bias in the parameter estimation 

Assuming the measurements of the galaxy power spectrum and bispectrum are Gaussian distributed, we can 
compute the likehhood of the measurements given a set of model parameters p, 

2£(p) =lndctC + (x-/x)'^C"^(x-/x) + constant, (B.l) 

where x and /i.(p) respectively represent the measurement and the model prediction of the galaxy power 
spectrum Pg(fc) and bispectrum _Bg(ka, k;,, kc), and C(p) is the covariance matrix for the measurements. For 
forecasting the parameter constraints in galaxy surveys in section 5 we have used the Fisher information matrix 
(see, e.g., [37, 38]) 

where the first term involving the derivative of the covariance matrix is smaller than the other, and hence we 
ignored the derivatives of the covariance matrix in computing the constraints. 

When the assumed fiducial model is incorrect, our computation of the mean and its covariance matrix 
will misrepresent the measurements, providing biased best-fit parameters p different from the true parameters 
■pt- Assuming that our model is not far off from the true model, the best-fit parameters can be approximated 
as p = p4 + (5p, and the bias (5p in the parameter estimation can be obtained from the relation that the best-fit 
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parameters p maximize the likelihood in eq. (B.l), 

' dL 







dpi 



(p) ) = Tr [C-^C., - C-iC.,C-i(D) + C-\T>, 



(B.3) 



Tr 



C-iC,,C-^(C - Cf) - 2J/x^C-V,t +2^(5pjTr -C-^C^.C-iC j + /x^^C" V,, 



The data matrix is D = (x — /x)(x — /x)"^, and the ensemble averages of the data matrix and its derivative are 

(D) = Ctivt) + (Mt - ^Ct-Y, C.,5p, + 0{Sp^) , (B.4) 

j 

(D,,) = ^ ti^^6p, - 2/x^,^/x^ , (B.5) 



where Ct is the true covariance matrix and the mean is related to the true mean as /x(pt) = /x(p) + Sfi. 
Substituting these relations in eq. (B.3), we can obtain the bias in our parameter estimation as 



ic-iC,,C-i(Q - C) + ^/x^C-V., 



^ cr| dpi ' I ^ 



a\ dpi dpj 
(B.6) 



where we again ignored the derivative of the covariance matrix. 
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